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Abstract 

The recently published experimental dependence of the J/tp suppression pat- 
tern in Pb+Pb collisions at the CERN SPS on the energy of zero degree 
calorimeter Ezdc are analyzed. It is found that the data obtained within 
the 'minimum bias' analysis (using 'theoretical Drell-Yan') are at variance 
with the previously published experimental dependence of the same quantity 
on the transversal energy of neutral hadrons Et- The discrepancy is related 
to the moderate centrality region: 100 < N p < 200 (N p is the number of 
nucleon participants). This could result from systematic experimental errors 
in the minimum bias sample. A possible source of the errors may be con- 
tamination of the minimum bias sample by off-target interactions. The data 
obtained within the standard analysis (using measured Drell-Yan multiplicity) 
are found to be much less sensitive to the contamination. 
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Recently, the NA50 collaboration published new data on the centrality dependence of 
the J/ip suppression pattern in Pb+Pb collisions at CERN SPS [[jj]. The centrality of the 
collisions was estimated by measuring the energy of projectile spectator nucleons E ZDC by 
a zero degree calorimeter. The purpose of the present comment is to compare the new data 
with the ones published by the NA50 collaboration previously where the transverse 

energy of produced neutral hadrons Et was used as the centrality estimator. It is shown 
that the two sets of data (for brevity, we shall call them '.E^c'-data' and 'E+-data', respec- 
tively), having very similar qualitative behavior, seem, nevertheless, to be at variance on 
the quantitative level. This may point out a presence of systematic errors in one of the data 
sets. 

In each data set, the J/ip suppression pattern were obtained in two different ways: by 
the standard analysis and the minimum bias analysis. In the standard analysis, the ratio 
R of J/ip multiplicity to the physical multiplicity of Drell-Yan pairs was calculated from 
the measured dimuon spectra. In the minimum bias analysis, the ratio R was obtained by 
dividing the experimentally measured ratio of the number of J /if) events to the number of 
minimum bias events by the theoretical Drell-Yan' given by 

mv(w u ^d%P(E T \b){DY{b))P int (b) 

{DY{Et)) = J* d%P(E T \b)P int (b) (1) 

Here (DY(b)) is the average number of Drell-Yan dimuon pairs per Pb+Pb collision at 
impact parameter b. Pi n t{b) is the probability that an interaction between two nuclei at 
impact parameter b takes place. Both quantities were calculated in the Glauber approach. 

The .E T -distribution of events at fixed impact parameter b, P(E T \b), was assumed to 
have the Gaussian form with the central value and dispersion given by 



(E T (b)) = qN p (b) a Er (b) = q^JaN p (b) (2) 

Here N p (b) is the average number of nucleon participants in a Pb+Pb collision at impact 
factor b calculated in the Glauber approach. The parameter values q = 0.274 GeV and 
a = 1.27 || were fixed from the minimum bias transverse energy distribution ||. 

The -EzDC-dependence of 'theoretical Drell-Yan' was found from the formula similar to 
(H). The .E'z.DC'-distribution P[E ZDC \b) was assumed to be a Gaussian with the central value 
and dispersion given by Eqs. (1) and (2) of Ref. 

It is obvious that in absence of essential systematic errors, the Et- and EzDC-d&t& should 
be consistent with each other. If there exists a parameterization of the 6-dependence of the 
J/ip multiplicity (J/ip(b)), such that the Et dependence of R computed from 

p ,p \ ^d^bPiET^iJ/m)^) 

r{Et) = Wh T (3) 

agrees with the .E^-data, then the dependence of R on E ZDC 

R(F S ^d%P{E ZDC \b){J/m)P l ntib) 

R{E ZDC ) = — — (4) 

\ UY )E ZDC 

should agree with the ^DC-data without any additional parameter fit. And vice versa: a 
model fitted to the _Ez,Dc-data should automatically agree with the E^-data. 
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To check, whether the two sets of the NA50 data indeed possess the above property, we 
used three different models to parameterize (J/if)(b)). Two of them, the models A and B (see 
below), were fitted to the E T -data presented in Fig. 4 of Ref. |4j] and then were compared 
with the 'minimum bias' EzDC-d&ta, from Fig.5 of Ref. Q. The third one, the model C, was 
fitted to the 'minimum bias' EzDC'-dsita. and compared with the E'r-data. 

Model A is essentially the statistical coalescence model (SCM) [§,0, which assumes 
that charmonia are created at the final stage of the heavy ion reaction from the charmed 
quark- ant iquark pairs produced in hard parton collisions at the initial stage. The free 
parameters of the SCM (see Ref. || for details) are a^ N , the in-medium modified || cross 
section of open charm production in nucleon-nucleon collisions, and rj, the fraction of dimuon 
pairs, originating from J /if) decays, that satisfies the kinematical conditions of the NA50 
spectrometer. 

The SCM is known to be applicable only for sufficiently large systems (N p > 100). To 
parameterize the data in the low E? region, we assumed a smooth transition between the 
normal nuclear suppression (NNS) to the statistical coalescence model: 

( JM&)> = exp [-(N 2 p (b)] (J/m)NNS + {l - exp [-(Nfa)] } (J/if)(b))scu. (5) 

This introduces an additional free parameter (. The value of NNS cross section was fixed 
at a a bs = 6.4 mb. 

The drop of R(Et) at Et ~ 100 GeV appears in SCM due to fluctuations of the ther- 
modynamical parameters and due to Et losses in the dimuon event sample relative to the 



minimum bias one 10 



A good fit of the Sy-data (x /dof = 1.10) is obtained with 

aj™ = 3.67 yb rj = 0.129 ( = 1.04 • 10~ 4 . (6) 
Model B is the geometrical model of charmonium suppression |ll| , [12| . It is assumed 



that all charmonia suffer the normal nuclear suppression with the cross section a a b s . Then 
all excited charmonium states are destroyed in the region, where the nucleon participant 
density in the plane transverse to the collision axis exceeds the value ri\. This suppresses 
the J /if) yield by a factor X < 1, because excited charmonia contribute to the final J /if) 
multiplicity through their decays. No charmonium state survives in the region, where the 
participant density exceeds the value Fluctuations of Et result in additional suppression 
at Et ^ 100 GeV W\. E T losses are ignored. 



To make our parameterization as good as possible we treated not only n\ and ri2, but 
also cr a b s and X as free parameters. A very good fit (x 2 /dof = 0.95) is obtained at 

m = 2.97 far 2 n 2 = 4.05 far 2 ( 
a abs = 5.16 mb X = 0.366 ^ ' 

In spite of the fact that the underlying ideas of the models A and B are very different, 
the fits are completely consistent with each other everywhere, except the low Et region, 
where only two points with large statistical errors are available (see Fig. 0). 

The both parameterizations are, however, in obvious disagreement with the 'minimum 
bias' i?2_DC*-data: x 2 /dof = 12.1 and x 2 /dof = 11.5 for the models A and B, respectively. 
The strongest discrepancy is observed in the domain 100 < N p < 200 (see Fig. §). 
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It must be noted that the effects that were responsible for the drop of the ratio R at 
E T > 100 GeV, the fluctuations and E T losses, do not influence the E Z dc dependence. 
Therefore, neither of the models is able to reproduce the drop of R(E Z dc) at E Z dc ~ 9 
TeV. 

Let us go the opposite way, i.e. we fit the Ezdc data and check, whether the fit agrees 
with the Et data. For this purpose, the model C will be used. 

Model C was proposed by NA50 collaboration in Ref. The idea is similar to the 
model B, but the parameter that controls the charmonium suppression is not the nucleon 
participant density, but rather the number of participants. It is assumed that at N p > N p i, 
the J/i/j yield is suppressed by the factor of X\ < 1 relative to the normal nuclear suppression 
value. At N p > N p2 the remaining J/ip yield is additionally suppressed by the factor X 2 . A 
minimum of x 2 is reached at 

N pl = 122 N p2 = 335 a abs = 5.47mb 

X 1 = 0.670 X 2 = 0.752 1 j 

Although the quality of the fit is not very good (x 2 /dof = 2.73), the model satisfactory 
reproduces the shape of the Ezdc data, including the drop of the ratio R at Ezdc ~ 9 TeV. 

Comparing the model C with the Et data reveals strong disagreement: x 2 /dof = 13.3. 
The strongest discrepancy is again related to the region 100 ^ N p ^ 200 (see Fig. [I]). At 
large E T , although the drop of R at E T ^ 100 GeV is present on the curve, it does not 
reproduce the behavior of the data. Account for Et losses in the dimuon sample (similar to 
[13|,10]) somewhat improves agreement with the data, but still the shape of the drop is not 



reproduced. 

The above consideration demonstrates that the two sets of the NA50 data on the cen- 
trality dependence of the J/ ip suppression pattern: the ^-dependence in Fig. 4 of Ref. (§] 
and the E^Dcrdependence in Fig. 5 of Ref. |jj are at variance with each other. This could 
mean that at least one of the two sets may contain sizable systematic errors. 

In our opinion, the E^DC-data obtained within the minimum bias analysis are very likely 
to contain essential systematic errors because of possible contamination of the minimum bias 
sample by off-target interactions. In the data analysis procedure of Ref. [TJ, empty target 
runs were used to remove the contribution of the off-target interactions. This approach, 
however, does not take into account possible influence of target on the shape of the con- 
tamination. In fact, spectator fragments from an off-target collision may interact with the 
target. On the other hand, the spectator fragment from a normal collision with a target 
nuclei may experience off-target interactions before they are registered by the zero degree 
calorimeter. Therefore, in presence of the target, the off-target contribution may be shifted 
downwards in E Z dc- 

Indeed, the minimum bias Et distribution perfectly agrees with the Glauber model (see 
Fig. 1 of Ref. In contrast, the agreement of the Ezdc minimum bias distribution seems 
to be not so perfect. Unfortunately, the authors of Ref. do not quote the value of x 2 /dof, 
which characterizes the quality of the fit of the minimum bias Ezdc distribution (shown in 



Disagreement is observed only in the low Et region, because the efficiency of the target algorithm 
is less than unity there. 
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Fig.l of Ref. [|TJ). Nevertheless, it is seen from the plot, that the experimental points lie 
slightly above the theoretical curves in the region 25 <^ Ezdc ~ 18 TeV (which corresponds 
to 100 < N p < 200). Although the difference does not exceeds even the size of the point 
symbols, it, nevertheless, indicates a sizable discrepancy, because of the logarithmic scale of 
the vertical axis of the plot. Note, that at low Ezdc, i>e. in the region of the 'second drop' 
of the Ez DC-data, the theoretical models are also in disagreement with the data. Therefore, 
the discrepancy of the two sets of data in the region 100 ^ N p ^ 200 and the drop of 
the ratio R(Ezdc) a ^ Ezdc ~ 9 TeV may have a common origin: contamination of the 
minimum bias sample by off-target interactions. 

The standard analysis data appear to be less sensitive to the off-target contamination. 
In fact, the models A and B are consistent with the standard analysis £xDC-data (presented 
in Fig. 3 of 0): x 2 /dof = 0.85 and x 2 /dof = 0.83, respectively (see Fig. |]). No definite 
conclusion is possible because of the large statistical errors. The model C may be also 
considered to be consistent with the standard analysis data with slightly worse but still 
quite satisfactory fit quality: x 2 /dof = 1.10. Nevertheless, no contradiction between E T and 
Ezdc is seen in the recently presented by the NA50 collaboration standard analysis data 
from the year 2000 run, which have essentially smaller error bars [fH ]. This allows to conclude 
that the minimum bias sample seem to be more sensitive to the off-target contaminations 
than the dimuon one. 

The discrepancy between E T and E Z dc data are unlikely to be related to interactions 
of Pb nuclei and their fragments with air. In fact, comparing of the data collected in the 



year 2000 with the target in vacuum |L4j with the data from previous runs |]-|J indicates 
that interactions with air influence both Et and Ezdc data and that the influence is related 
mostly to the peripheral domain N p ^ 100. 

The contaminations at moderate centrality, 100 < N p < 200, may come from interactions 
of reactions fragments with nuclei in zero degree calorimeter. If such an interaction takes 
place, a fraction of the produced particles (those having large transversal momenta) may 
leave the calorimeter without deposing their energy in it. This could cause sizable losses 
of zero-degree energy, which result in distortion of the shape of the minimum bias sample. 
Therefore, the mentioned problem will likely persist even in the new vacuum data, if they 
are analyzed within the 'minimum bias' procedure, unless losses of the zero degree energy 
are taken into account. 

We conclude, that the new data on the Ezdc dependence of the J/xjj suppression pattern 
in Pb+Pb at the CERN SPS [|I|] are at variance with the other data published by the same 
collaboration [fj-|4j]. This problem is related only to the data obtained within the 'minimum 
bias' analysis. A possible source of the discrepancy might be distortion of the minimum 
bias sample most likely by losses of zero degree energy because of interactions of reaction 
fragments with nuclei inside the zero degree calorimeter. 
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FIG. 1. The dependence of the J/ijj to Drell-Yan ratio R on the transverse energy Et- The 
corresponding average number of participant nucleons N p are shown on the upper axis. The points 
with error bars are the NA50 data. The lines correspond to different parametrizations. Models 
A an B are fitted to the Et data set. Model C being fitted to the Ezdc data set (see Fig. ||) 
demonstrates disagreement with the Et data. 
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FIG. 2. The dependence of the J/ip to Drell-Yan ratio R on the energy of zero degree calorimeter 
Ezdc- The corresponding average number of participant nucleons N p are shown on the upper axis. 
The points with error bars are the NA50 'minimum bias' analysis data. The lines correspond to 
different parameterizations. Model C is fitted to the Ezdc data set. Models A and B being fitted 
to the -Er-data (see Fig. ||) demonstrate disagreement with the Ezdc data. 
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FIG. 3. The dependence of the J / to Drell-Yan ratio R on the energy of zero degree calorimeter 
Ezdc- The corresponding average number of participant nucleons N p are shown on the upper axis. 
The points with error bars are the NA50 standard analysis data. The lines correspond to different 
parametrizations. Models A and B are fitted to the -Er-data (see Fig. [I]) Model C is fitted to the 
'minimum bias' Ezdc data (see Fig. |2|). Due to rather large statistical errors, all three models 
may be considered to be consistent with the standard analysis data. 
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